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Abstract 

Results are presented from a phenomenological analysis of recent measurements of jet suppression and 
modifications of jet fragmentation functions in Pb+Pb collisions at the LHC. Particular emphasis is placed 
on the impact of the differences between quark and gluon jet quenching on the transverse momentum (pC) 
dependence of the jet Raa and on the fragmentation functions, D{z). Primordial quark and gluon parton 
distributions were obtained from PYTHIAS and were parameterized using simple power-law functions and 
extensions to the power-law function which were found to better describe the PYTHIAS parton spectra. 
A simple model for the quark energy loss based on the shift formalism is used to model Ran and D{z) using 
both analytic results and using direct Monte-Carlo sampling of the PYTHIA parton spectra. The model is 
capable of describing the full rapidity, and centrality dependence of the measured jet Raa using three 
effective parameters. A key result from the analysis is that the D{z) modifications observed in the data, 
excluding the enhancement at low-x, may result primarily from the different quenching of the quarks and 
gluons. The model is also capable of reproducing the charged hadron Raa at high transverse momentum. 
Predictions are made for the jet Raa at large rapidities where it has not yet been measured and for the 
rapidity dependence of D{z). 


1. Introduction 

Measurements of jet production and jet proper¬ 
ties in ultra-relativistic nuclear collisions provide 
an important tool to study the properties of quark 
gluon plasma created in the collisions. High-energy 
quarks and gluons produced in hard-scattering pro¬ 
cesses can interact with and lose energy while prop¬ 
agating in the plasma. Those interactions can 
both reduce the energy of the jets that result from 
the fragmentation of the quarks and gluons and 
change the properties of the jets. These and other 
“medium” modifications of the parton showers ini¬ 
tiated by the hard scattering [1, 2] are frequently 
collectively referred to as “jet quenching”. 

Jet quenching was first observed at the LHC 
through the observation of highly asymmetric dijet 
pairs [3] that result when the two jets lose different 
amounts of energy in plasma. Since dijet pairs for 
which both jets lose similar energy or, more gener¬ 
ally, have similar modifications will appear “sym¬ 
metric” , other observables are needed to probe the 
effects of quenching on the typical jet. Measure¬ 
ments of the suppression of the hadron spectrum 
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resulting from the energy loss of the parent jets 
have been carried out at both RHIC [4-6] and the 
LHC [7-9]. These show a suppression that at the 
LHC varies from a factor ^ 5 for hadron trans¬ 
verse momentum (px) values ^10 GeV to a fac¬ 
tor of ~ 2 for pt L 50 GeV. Most jet quenching 
calculations that attempt to infer medium proper¬ 
ties such as the quenching transport parameter, q, 
(see e.g. [10] and references therein) have relied on 
the single hadron suppression results because of the 
theoretical simplicity in calculating single hadron 
spectra. However, the single hadron measurements 
have only indirect sensitivity to the kinematics of 
the parent parton and little sensitivity to the details 
of the modification of the parton shower. 

Recent measurements of the suppression of the 
jet yield in 2.76 TeV Pb+Pb collisions at the LHC 
[11] are expected to provide a more sensitive probe 
of the physics of jet quenching at least through 
the improved correlation between the measured jet 
and the parent parton (shower) kinematics. Recent 
measurements of the jet nuclear modification fac¬ 
tor, Raa, for high transverse momentum jets show 
a factor of ~ 2 suppression in the jet yield that 
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increases slowly with increasing The suppres¬ 
sion is observed to vary monotonically as a func¬ 
tion of collision centrality and to be independent 
of jet rapidity within the statistical and system¬ 
atic uncertainties. Separately, recent measurements 
of the fragmentation functions for jets produced in 
Pb-|-Pb collisions [12, 13] have shown an enhanced 
yield of hadrons with low transverse momenta, pt, 
or low jet momentum fraction, z, a suppressed yield 
of hadrons with 0.05 < z < 0.2 and possibly an en¬ 
hanced yield of hadrons with z > 0.2. These two 
observables, the single jet suppression and the mod¬ 
ification of the fragmentation function, arguably 
provide the minimal set of data needed to under¬ 
stand the physics of jet quenching: the jet suppres¬ 
sion is sensitive to the amount of energy the jet loses 
- outside the jet “cone” - while the fragmentation 
function is sensitive to the re-distribution of energy 
inside the jet cone. 

Interpretation of the inclusive jet suppression and 
fragmentation data is complicated by the flavor ad¬ 
mixture of the primordial partons. For this anal¬ 
ysis, we focus on the relative combination of light 
quarks and gluons which are expected to suffer dif¬ 
ferent energy loss due to their different color charges 
and/or differences between the quark and gluon 
splitting functions. In weak coupling calculations, 
the relative quark and gluon energy loss rates are 
determined by perturbative QCD color factors. For 
example, it is usually assumed that gluons lose en¬ 
ergy at a rate 9/4 higher than that for quarks. 
Recent studies including those based on quench¬ 
ing Monte Carlo codes have been compared to the 
jet i?AA and fragmentation function measurements, 
but these analyses have not explicitly attempted to 
elucidate the role of the relative quark and gluon 
contributions to the jet spectrum. 

In the remainder of this paper, we attempt to 


interpret recent measurements of the single jet sup¬ 
pression and fragmentation function ratios explic¬ 
itly accounting for the role of the quark and gluon 
admixture. We argue that some of the features of 
the data can be explained purely on the basis of 
the pi^* dependence of the quark to gluon fraction 
of the primordial parton spectrum and the different 
quark and gluon energy loss. Our analysis is based 
on simple assumptions regarding the parametric de¬ 
pendence of the energy loss on the jet transverse 
momentum and flavor. These assumptions are suf¬ 
ficiently simple that the results of our analysis can 
be easily explained and understood, but their sim¬ 
plicity also means that the results presented here 
should be verified using a proper jet quenching cal¬ 
culation. Our analysis is based on a combination 
of analytic calculation and Monte Carlo simulation 
using simulated quark and gluon spectra obtained 
from PYTHIAS [14]. 

2. Parameterizing Jet Spectra and D(z) dis¬ 
tributions 

A key ingredient of the analysis in this paper 
is the quark and gluon jet spectra which were ob¬ 
tained from PYTHIAS using procedures chosen to 
be similar to those used in certain ATLAS simu¬ 
lations [15, 16]. Namely, PYTHIAS was run using 
parameters from the AU2 tune [17] and using CTIO 
parton distribution functions [IS]. This combina¬ 
tion was shown to describe well the LHC jet data 
[16]. Jets were reconstructed using the anti-fcj al¬ 
gorithm [19] with the distance parameter R = 0.4 
applied to hadrons with lifetimes cr > 1 mm. The 
resulting jets were matched to one of the two outgo¬ 
ing partons from the leading order hard-scattering 
process by choosing the parton with the smallest 
angular distance, AR = ^+ Arj'^ and the fla- 


Fit type 

Parameter 

|y| <2.1 

\y\<0.3 

0.3 < \y\ < O.S 

1.2< [yj <2.1 

All 

fia 

0.34 

0.2S 

0.29 

0.40 

Power law 

riq 

5.66 

5.37 

5.40 

6.15 

Ug 

6.25 

5.97 

6.09 

6.92 

Extended 
power law 

Uq 

4.19 

4.34 

4.27 

3.75 

P, 

0.71 

0.49 

0.54 

1.2 

Ug 

4.69 

4.55 

4.57 

4.60 

Po 

O.SO 

0.71 

0.76 

1.2 


Table 1: Parameters obtained from fits of the PYTHIAS jet spectra to power-law (Eq. 1) and extended power-law (Eq. 2) 
functions. 
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Figure 1; Top panels: Quark (left) and gluon (right) spectra obtained from PYTHIAS simulations over the rapidity interval 
\y\ < 2.1 with simple power law (Eq. 1) and modified power law (Eq. 2) fits superimposed. Bottom panels: ratios of the spectra 
to the two fit functions. 


vor of the jet was assigned to be that of the matched 
parton. Six million PYTHIA hard-scattering events 
were generated for each of five intervals of ptj the 
transverse momentum of outgoing partons in the 
2 —I 2 hard-scattering, with boundaries 17, 35, 70, 
140, 280 and 560 GeV. 

The Monte Carlo results presented in this paper 
were obtained by directly using the jets obtained 
from the PYTHIAS simulations, but the analytic 
results require a parameterization of the depen¬ 
dence of the jet yields. A common parameterization 
used to describe the spectra of high-px hadrons or 
jets is the power-law form, e.g. 


dn 



( 1 ) 


where pxo is a reference transverse momentum 
value at which A represents dn/dp^^*". Results of 
fits to the quark and gluons distributions from the 
IpI < 2.1 rapidity interval are shown in Fig. 1. The 
values of n extracted from the pure power-law fits 


are listed in Table 1. The power-law function can 
describe the gross-features of the spectra but the 
ratios of the spectra to the fit functions presented 
in the bottom of the figure indicate significant de¬ 
viations of the jet spectra from the power-law form. 

The power-law distribution can be improved by 
adding a logarithmic pl^* dependence to the expo¬ 
nent producing an “extended power-law”, 

V n+/31og (p^^Vpto) 

j ■ (2) 

With this form, j3 represents the logarithmic deriva¬ 
tive of dn/dp!^^ at pl^* = pxq- At the most forward 
rapidities, the strong phase-space suppression of the 
jet spectra at high pl^* makes even the extended 
power-law inadequate for describing the jet spec¬ 
tra. Thus, for the most forward rapidities, an ad¬ 
ditional quadratic term, qlog^ (p^t^/pto^ > added 
to the power-law exponent and the resulting func¬ 
tion is capable of describing the most forward quark 


djff ( 
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and gluon spectra over the range used in this 
analysis. 

A jet spectrum that consists of a mixture of quark 
and gluon contributions can be represented in terms 
of a sum of contributions each of the form of Eq. 1 
or its extensions. However, for the purposes of this 
paper, it will be convenient to express the combined 
spectrum in terms of a quark fraction, /^p, specified 
at pto- Then a combined spectrum using power-law 
forms can be written 


( 3 ) 

where and ng are the quark and gluon power- 
law indices, respectively. Since Ug ^ rig^ the quark 
fraction will evolve as a function of pl^* according 
to 


dN 


= A 




Pto 



'90 




'90 



+ (1 “ flo) 


1 



1-h 


1-/90 

/90 



( 4 ) 


For the extended power-law parameterizations of 
the spectra, the pl^'^-dependent quark fraction looks 
similar to that in Eq. 4 but with the addition of a 
term, (/3g —/3g)log ^pl^^/pTo'j to the exponent in 

the denominator. The pl^* dependence of the quark 
fraction is shown in Fig. 2. 

The PYTHIAS D{z) distributions were obtained 
using final-state charged hadrons located within an 
angular radius, AR < 0.4, of reconstructed jets 
having pl^* > 100 GeV. The resulting distribu¬ 
tions are shown in Fig. 3 for the rapidity interval 
\y\ < 2.1. The quark D{z) distribution is notice¬ 
ably harder than the gluon D{z) distribution, but is 
also lower at intermediate z, in the range where the 
D{z) distribution appears to be depleted in Pb-|-Pb 
collisions. 

For use in the analytic analysis, the D{z) distri¬ 
butions were fit to functions of the form, 

“ ■ 1! T 

(1 -I- ezj'^ 

which are similar to other commonly used parame¬ 
terizations [20] with the addition of an exponential 



different jet rapidity intervals used in this study. The points 
show results obtained from PYTHIAS simulations, the solid 
lines represent results obtained from extended power-law fits 
with the parameters shown in Table 1. 



a 

h 

c 

d 

e 

/ 

quark 

31S 

2.51 

1.44 

-0.S5 

52.4 

0 

gluon 

574 

1.S7 

2.32 

9.09 

32.0 

10.3 


Table 2; Parameters describing the fragmentation functions 
extracted from PYTHIAS using the procedures described in 
the text for the functional form in Eq. 5. 


term. That term is not used for the quark distribu¬ 
tions, but it’s presence provides a more controlled 
description of the gluon D{z) distribution. The re¬ 
sults of the fits for the quark and gluon distributions 
over \y\ < 2.1 are shown in Fig. 3, and the ratios 
of the fit to the PYTHIAS D{z) distributions are 
shown in the lower panels. The fits well describe the 
simulated D{z) distributions with parameters that 
are provided in Table 2. We note that the parame¬ 
terization in Eq. 5 has a smooth extrapolation past 
z = 1. The pQCD fragmentation function has no 
contribution from z > 1, but when reconstructing 
jets in PYTHIAS and data, there are events having 
two jets that are close enough that a high-plj/ frag¬ 
ment from the higher-energy jet can be associated 
with the lower-energy jet possibly yielding a hadron 
with z > 1. The D{z) distributions fall rapidly 
above z = 1 so they have no practical importance, 
though the continuity of the parameterization will 
be relevant later in this paper. 
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the analog of the measured Raa, is then 



Figure 3: PYTHIAS quark and gluon D(z) distributions 
for R = 0.4 jets having > 100 GeV and \y\ < 2.1. The 
solid lines show the results of fits to the D{z) distributions 
using the function in Eq. 5. 

3. Analytic models for jet Raa and D{z): 
power-law spectra 

Our modeling of the single jet suppression is 
based on the “shift” approach of [21] in which 
the quenched jet spectrum can be (approximately) 
written 

dpi^* y dPr*/ 

( 6 ) 

where dnq and dn represent the per-event yields of 
quenched and unquenched jets. The second term in 
Eq. 6 is a Jacobian term that is necessary to (e.g.) 
preserve the total number of jets. 

Using the power-law form for a single-flavor jet 
spectrum, 

= ... + 
dp^T XP^T+^ip'r)) \ dp!^ j 

(7) 

The ratio of the quenched and unquenched spectra, 


7?aa(p^t*) 


(l + 5(pir^‘)/pir^‘ 


n 



dS \ 


( 8 ) 


It has been previously observed [22] ^ that if the 
shift is proportional to pl^*, S = spx, then the re¬ 
sulting Raa is pl^*-independent: 


7?aa(p¥‘) = ... , \n-i ^ (9) 

[l + s) 


such that the fractional shift can be inferred from 
an approximately pl^*-independent Raa [22] 


1 


( 10 ) 


This result has previously been applied to single 
hadron Raa measurements but is arguably more 
appropriate when applied to jet Raa measure¬ 
ments. Then, naively applying Eq. 10 to the typical 
suppression observed at high pl^* in central Pb-|-Pb 
collisions at the LHC, Raa ~ 0.5, and using n ^ 5, 
the resulting fractional shift would be s = 0.19. 
The observation that the jet Raa is only weakly 
dependent on pl^* at high pl^* has been taken as 
evidence that jets lose a constant fraction of their 
energy in the quark gluon plasma created in Pb-l-Pb 
collisions. There are potentially significant theoret¬ 
ical flaws with this conclusion, but the conclusion 
and the interpretation of the extracted s also suf¬ 
fer from neglecting the fact that the jet spectrum 
is composed of an admixture of flavors. 

Starting with a combination of quark and gluon 
power-law spectra (Eq. 3), the quenched spectrum 
would be 


dp!^ 


Pto 


igo 


.p't 



dsA 



rig 



( 11 ) 


The resulting Raa would be given by the ratio of 
Eq. 11 to Eq. 3 which, with some simplification, 


^The definition of n differs between this paper and [22] 
where it characterizes the invariant cross-section and is, 
therefore, larger by one due to the dp^ factor. 
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takes the form 


1 


loss. Starting from Eq. 11, the modified quark frac¬ 
tion in the constant fractional shift scenario is 


Raa = fq 


1 + Sq/p!^ 


jet 




dd 


jet 


(1 - fq) 


dSa 


, 1 + 


jet 


djf^ 


jet 


( 12 ) 


Here, /q is the full pl^*-dependent quark fraction 
in Eq. 4. As the equation indicates, the combined 
jet Raa is given by a combination of the separate 
quark and gluon suppression factors weighted by 
the quark and gluon fractions. 

Even in the case of constant fractional energy 
loss for the quarks and gluons, Sq = Sqp^^^ and 
Sg = 12 would imply an Raa that varies 

with as long as Sq ^ Sg. Eor example, if 
Sg/Sq > 1, then the Raa will increase with 
because of an increasing quark fraction and weaker 
suppression for the quarks. Such an increase of Raa 
with increasing pl^* has been observed in the mea¬ 
sured Raa values that are shown in Fig. 4. To test 
whether the data are compatible with the constant 
fractional energy loss scenario, we have assumed 
Sg = 9/4 X Sg and have fit the Raa values using 
Eq. 12 with one free parameter for each centrality 
bin, namely Sq. The results are shown in Fig. 4 with 
solid lines. The extracted Sg values vary from 0.02 
for the 60-80% centrality bin to 0.1 for the 0-1% bin. 
The figure shows that the constant fractional shift 
assumption combined with the power-law form for 
the jet spectra is capable of approximately repro¬ 
ducing the slow variation of the measured jet Raa 
with pi^*. 

Since the pl^* dependence of the jet suppression 
can be successfully explained using a combination of 
the varying quark fraction and the greater quench¬ 
ing of gluon jets, it is worth exploring the impact of 
these same behaviors on the jet fragmentation func¬ 
tion. In particular, we wish to determine whether 
the different quenching of the quarks and gluons 
can explain part or all of the the observed modifi¬ 
cations of the fragmentation function in Pb-|-Pb col¬ 
lisions. We make the simplest possible assumption, 
namely that the quarks and gluons lose energy in 
the plasma and then fragment according to vacuum 
fragmentation functions. With this assumption, the 
only source of modification to the inclusive jet frag¬ 
mentation function is the change in the quark (or 
gluon) fraction due to the medium-induced energy 


fq 


mod 


1-^ 



(1 + ^.)"-'^ 

{l + Sq^-^ Ip¥V 


(13) 

Assuming that the D{z) distributions are inde¬ 
pendent of pl^*, the per-jet distribution of fragments 
as a function of the fragment longitudinal momen¬ 
tum fraction, D{z), can be written. 


D{z) = f^^Dqiz) + (1 - /f)i?g(z), (14) 


where Dq{z) and Dg{z) are the quark and gluon 
D{z) distributions, respectively, and /“* is the 
modified quark fraction integrated over a given pl^* 
range. 

The ATLAS jet fragmentation measurements 
were obtained for pl^* > 100 GeV. Applying Eq. 14 
over this pl^* range and using the Sg parameters 
obtained from hts to the jet i?AA, we calculated 
the ratio of modified D{z) distributions in differ¬ 
ent centrality bins to the distribution in the 60- 
80% centrality bin for comparison with the ATLAS 
data. The results are shown along with the data in 
Fig. 5. The hgure shows that our simple model for 
the medium modihcations of the inclusive jet frag¬ 
mentation function can reproduce some of the qual¬ 
itative features in the data, namely the suppression 
of the fragmentation function at intermediate z and 
an enhancement in the fragmentation function at 
large z. This latter is statistically marginal in the 
data given the (combined) error bars, but the en¬ 
hancement at large z in the model is an automatic 
result of the increased quark content of the jet spec¬ 
trum. Our model does not show as deep a suppres¬ 
sion in the D{z) ratio near 0.1 which may indicate 
that additional physics contributes there. 

One feature in the data that cannot be explained 
by the model is the enhancement at low z. Our sim¬ 
ple model also explains the centrality dependence of 
the data, except for the 50-60% centrality bin, given 
the fits to the single-jet suppression. Based on the 
results shown in Fig. 5 we argue that it is plausible 
that the modifications observed at intermediate and 
large z in the jet fragmentation function result from 
quenching-driven changes in the jet quark fraction 
while the enhancement at low z reflects a contribu¬ 
tion of extra particles in the jet either from radiative 
emission within the jet or recoil of particles in the 
medium. 
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Figure 4: Nuclear modification factor of jets, -RaAi measured by ATLAS [11] (black markers) in four different centrality bins 
(rows) and four different rapidity regions (columns) compared to the analytic calculation (red line) and MC calculation (blue 
histogram) of the Raa using constant fractional shift and power-law spectra. 


We have performed a separate Monte-Carlo eval¬ 
uation of the single-jet suppression to check and 
improve on the results of the above analytic calcula¬ 
tions which are necessarily limited by assumptions 
regarding the shapes of the jet spectra. To simu¬ 
late the single-jet suppression, we sample jets from 
the PYTHIAS-simulated events, apply the shift as 
in Eq. 6 with chosen Sg and Sg for quark and 
gluon jets, respectively, and then build the result¬ 
ing spectra of quenched jets. The simulated Raa 
is obtained from the ratio of the quenched spec¬ 
trum to the original spectrum of PYTHIAS jets. 
The results are shown with the blue histograms in 
Fig. 4. The agreement with the analytic results is 
poor, suggesting that the power-law parameteriza¬ 
tion of the jet spectra is inadequate for the sim¬ 


ulation of the single-jet suppression. In fact, the 
Monte-Carlo sampled Raa decreases with increas¬ 
ing in contradiction with the general result that 
constant fractional shift combined with the increas¬ 
ing quark fraction should produce an Raa that in¬ 
creases with The decrease of the Monte-Carlo 
sampled Raa must necessarily result from the jet 
spectra being steeper, or having greater curvature 
than the power-law function or, equivalently, from 
the differences between the power-law fit functions 
and the PYTHIAS spectra seen in Fig. 1. 
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Figure 5: Ratios of D{z) distributions for six bins in collision centrality to those in peripheral (60-80%) collisions, 

F^('2:)|cent/-C^('2:)|60—80> measured by ATLAS for R = 0.4 jets [12] (black markers) are compared to the analytic calculation 
(red line) and MC calculation (blue histogram) of the same quantity in the fractional energy loss model. The analytic calcula¬ 
tion uses the power law parameterization of jet px spectra. 


4, Analytic models for jet Raa- extended 
power-law spectra 


To test this conclusion, we have extended the an¬ 
alytic analysis from the previous section to the case 
of extended power-law spectra. The analog of Eq. 7 
is 


dp!^* 


PTq 


n+/3 log [(Pt‘+'S)/pto 


1 + 


and the Baa is: 


dS 

dp^r^* 

(15) ' 


/ \ / \ 2/31og(l + S/pjr'*) 

V n+^log(l+S/p^^‘) 




(16) 


The logarithmic term in the exponent of the ex¬ 
tended power-law function not only produces a cor¬ 
responding logarithmic term in the exponent of 


1/(1 -I- S'/Px*) but it also generates an explicit de¬ 
pendence on pto/Pt*- Thus, even for a constant 
fractional shift, the Baa decreases with increasing 
for positive /3: 


/ \ 2/31og(l-|-s) 


J®* / 


1 -I- s 


1-1-1-/3 log (1-l-s) 


(17) 

For a spectrum consisting of both quarks and glu¬ 
ons, the Baa is given by an expression similar to 
Eq. 12 containing an fg- and (1 —/q)-weighted com¬ 
binations of Eq. 15 with different values for n and 
/3 for quarks and gluons. Given the fg distributions 
in Fig. 2, a constant fraction shift is unable to re¬ 
produce the measured increase of the jet Baa with 
increasing pl^*. In fact, a calculation of the Baa 
using the extended power-law form well reproduces 
the results of the Monte Carlo evaluation shown in 
Fig. 4. Thus, we conclude that the apparent suc¬ 
cess of the constant fractional shift scenario in Fig. 4 
using the analytic analysis is false and results from 
neglecting the deviations of the jet spectra from the 
pure power-law form. 



















































5. Modeling jet i?AAand £>( 2 ;): non-constant 
fractional shift 


The inability of the constant fractional shift as¬ 
sumption to explain the dependence of the mea¬ 
sured Raa suggests that the shift has to vary with 
more slowly than linearly. In fact, theoretical 
analyses of medium-induced energy loss would not 
be compatible with a constant fractional shift sce¬ 
nario. To proceed, we make the minimal extension 
of the analysis above and assume that the shift is 
proportional to an undetermined power of 


S = s' 



(18) 


Lacking any knowledge of the appropriate scale for 
the term in the parenthesis, we use the same refer¬ 
ence scale, ptqi used for parameterizing the spectra. 
Then, s', which has dimensions of energy or trans¬ 
verse momentum, represents the shift in transverse 
momentum for jets having pl^* = ptq- The result¬ 
ing Raa for a single jet spectrum and for combined 
quark and gluon spectra can be obtained using the 
procedures described above, in particular a com¬ 
bination of Eq. 15 weighted by the pl^*-dependent 
quark and gluon fractions; the formulas are not 
shown here for sake of brevity. 

The shift expression in Eq. 18 and the result¬ 
ing Raa were used to perform a fit to the ATLAS 
data to extract s' and a in different centrality bins. 
The fits were performed using the statistical un¬ 
certainties in the The results are presented in 
Eig. 6 where the top panels show contours in s'-a 
space for the 10-20% (left) and 70-80% (right) cen¬ 
trality bins. As the contours demonstrate, there is 
a strong correlation between the parameters which 
causes the obtained optimal a and s' values, shown 
in the lower panels of the figure as a function of 
A^part, to fluctuate. The a values, however, clus¬ 
ter around an average value of 0.55. To reduce the 
point-to-point scatter in the obtained parameters 
a was fixed to the value 0.55 and the fits were run 
again to extract s'. The values are shown in the 
lower right panel with the red points. The white 
and black circles shown on the y^ contour plots in¬ 
dicate, respectively the results of the free fits and 
the fits with a = 0.55. 

The s' values obtained using hxed a show an ap¬ 
proximately linear dependence on A^part and vary 
from ^ 1 GeV in the most peripheral bin to 
~ 5.5 GeV in the most central (0-1%) bin. The 


fact that s' extrapolates to a non-zero value for 
Vpart 0 may indicate that there is an additional 
contribution to the measured single-jet suppression 
present in even peripheral collisions. In fact, the 
free fits suggest a systematic rise in a for the two 
most peripheral bins which may arise from the same 
underlying physics. 

The Raa calculated using the results of the fit¬ 
ting procedure are shown in Eig. 7. Using the 
extended power-law parameterization of the quark 
and gluon spectra, the analytic and Monte Garlo 
results are in good agreement. The growth of the 
Raa with results from the fact that the shift 
increases with more slowly than linearly. Thus, 
the fractional energy loss decreases with increas¬ 
ing The agreement with the data is largely 
by construction since the parameters of the energy 
loss were obtained from the above-described fitting 
procedure. Nonetheless, our model is capable of de¬ 
scribing the available data with a single, centrality- 
independent value for a and a proportionality con¬ 
stant, s' that varies approximately linearly with 

Vpart ■ 

The rapidity dependence, or lack thereof, in the 
Raa arises from a cancellation between the rapid¬ 
ity dependence of the quark fraction, which in¬ 
creases with increasing rapidity (see Fig. 2), and 
the shapes of the quark and gluon spectra which 
become steeper with increasing rapidity. This can¬ 
cellation is illustrated in Fig. 8 which shows the sup¬ 
pression for the quark, gluon, and combined spectra 
in the \y\ < 0.3 (left) and 1.2 < |y| < 2.1 (right) ra¬ 
pidity bins. The difference between the quark and 
gluon suppression is greater in the 1.2 < \y\ < 2.1 
than in the |p| < 0.3 bin. The pl^ dependence is 
also much flatter in the higher rapidity bin. Yet, 
the combined suppression taking into account the 
pl^* dependence of /q in Fig. 2 is nearly the same 
in the two rapidity bins. 

The D{z) distributions calculated using the ex¬ 
tended power-law functions and the shift in Eq. 18 
are shown in Fig. 9. As with the Raa, the agree¬ 
ment between the analytic calculation and the 
Monte-Garlo sampled result is much better using 
the extended power-law descriptions of the primor¬ 
dial parton spectra. However, the D(z) modih- 
cations in the model are largely the same using 
the fractional and non-fractional shift parameter- 
izations. This lack of sensitivity to S'(pl^*) arises 
because the D{z) measurements are dominated by 
contributions from jets with pl^* ^ 100 GeV and 
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Figure 6: Top: x^/DOF as a function of a and s’ for the 10-20% (left) and 70-80% (right) centrality bins. The positions of 
the minima with a free (fixed) are indicated by the white (black) circles. Bottom: Parameters of non-fractional shift (Eq. 18) 
model, a (left) and s' (right), as a function of A^part obtained from fits of the resulting calculated Raa to the ATLAS data. 
The blue points indicate results for which both a and s’ are free parameters while the red points indicate the results of fits 
with a fixed to 0.55 (see text) shown on the left panel by the red line. The line on the right panel shows the result of a linear 
fit to s'(fVpart). 


because the D(z) modifications in the model pri¬ 
marily result from the difference between quark and 
gluon quenching for jets with similar transverse mo¬ 
menta. Thus, as long as the model reproduces the 
Raa near 100 GeV the D{z) modifications will be 
insensitive to the p!^ dependence of S. 

6. Rapidity dependence of the suppression 

The fraction of jets initiated by light quarks 
evolves as a function of the rapidity such that the 
probability that the jet is initiated by a quark is 


increasing with increasing rapidity. The steepness 
of the jet pt spectrum also evolves as a faction of 
the rapidity such that the px spectra of forward jets 
are steeper than the spectra of jets produced in the 
central region. Both of these features are demon¬ 
strated in Fig. 2 and in Table 1 of Sec. 2. Both fea¬ 
tures also influence the jet Raa, though they act in 
opposite directions. Nonetheless, it can reasonably 
be expected that the jet Raa will exhibit a differ¬ 
ent behavior in the forward region compared to the 
central region, or, equivalently, that the Raa will 
vary with rapidity at sufficiently large values. Thus, 
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Figure 7: Nuclear modification factor of jets, TJaAi measured by ATLAS [11] (black markers) in four different centrality bins 
(rows) and four different rapidity regions (columns) is compared to the analytic calculation (red line) and MC calculation (blue 
histogram) of the same quantity in the non-constant fractional energy loss model. The analytic calculation uses the extended 
power law parameterization of the jet px spectra that includes the logarithmic dependence of the exponent on jet px- 


it is clearly of interest to test the model presented 
in this paper by predicting the jet Raa in the for¬ 
ward region where it has not yet been measured. To 
do that, the jet Raa was calculated using the ana¬ 
lytic model in two bins of jet rapidity corresponding 
to those used by ATLAS or CMS [23, 24], namely 
2.1 < \y\ < 2.8 and 2.8 < \y\ < 3.5. In the later 
rapidity region, the jet pr spectra decrease approx¬ 
imately by four orders of magnitude in the region 
of jet pt between 40 GeV and 100 GeV. This steep 
fall-off of the spectra was found to be insufficiently 
described by the modified power-law, Eq. 2. To im¬ 
prove the parameterization, an additional quadratic 
term was introduced leading to the parameteriza¬ 


tion, 

/ \ u+/31og(pTVpTo)-l-7log^ (ptVpto) 

Ut / 

(19) 

which was found to describe the PYTHIA jet px 
spectra at the level of accuracy better then 10%. 
The resulting parameters and the quark fractions 
for the jet px spectra selected in the two rapidity 
regions are summarized in Table 3. 

The resulting analytic Raa was calculated using 
an extension of Eq. 16 to account for the quadratic 
term, and using the results from Sec. 3, namely 
a shift of the form of Eq. 18 with a = 0.55 and 
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Figure 8: Quark, gluon, and combined Raa vs for the 
\y\ < 0.3 (left) and 1.2 < \y\ < 2.12 (right) rapidity bins. 


Parameter 

2.1 < \y\ < 2.8 

2.8 < \y\ < 3.5 

riq 

5.5 

6.7 

Ug 

6.3 

7.4 

P, 

0.34 

-0.46 

Po 

0.52 

-1.19 

Iq 

1.2 

2.6 


1.5 

2.4 

fq 

0.60 

0.76 


Table 3; Parameters obtained from fits of the PYTHIAS for¬ 
ward jet spectra to the extended power-law (Eq. 19) forms. 

s^(-^part) as shown in Fig. 6. 

The predicted forward Raa is shown as a func¬ 
tion of in Fig. 10. A clear change in the trend 
of the Raa evolution with jet px can be seen. In 
contrast to the slow increase seen for the jet Raa 
in the rapidity regions within \y\ < 2.1, the jet 
Raa in the forward regions first increases, reaches a 
maximum and then decreases with increasing pl^*. 
The decrease is more pronounced for more forward 
region where the jet Raa in 0-10% central colli¬ 
sions reaches the maximum of approximately 0.4 
at around 50 GeV, and then it decreases reaching 
a value of approximately 0.15 at 170 GeV. These 
trends are present across different centralities. Such 
pronounced change in the behavior of the forward 
jet Raa represents a distinct feature that can be 
tested by future measurements at the LHC. 

The dependence of the quark fraction on the jet 
rapidity has to influence also the trends measured 
in the centrality dependent ratios of fragmentation 
functions, Rd(z), presented in Sec. 5. As demon¬ 
strated in Sec. 3, this ratio exhibits only a weak 


dependence on the shape of the underlying jet px 
spectra and thus, it is a very useful observable that 
may help isolate the effects of different quenching of 
quark-initiated and gluon-initiated jets. The sen¬ 
sitivity of the D{z) modification to differences in 
quark and gluon quenching are illustrated in the left 
panel of Fig. 11 which shows Rd{z) in 0-10% cen¬ 
tral collisions evaluated for three different choices 
of an effective color factor, Cp = 1.0,9/4, and 
3.0. A clear dependence of the Rd{z) on the color 
factor can be seen. For equal suppression of the 
quark-initiated and gluon-initiated jets, the Rd(z) 
exhibits only negligible difference from unity re¬ 
flecting minor difference in the quark fraction at 
unquenched and quenched jet transverse momenta. 
For color factor larger than the default value of 9/4, 
the Rd(z) exhibits larger increase compared to the 
Rd{z) evaluated with the default value of Cp. 

The middle panel of Fig. 11 shows the Rd{z) 
evaluated in 0-10% central collisions for the ra¬ 
pidity region |p| < 2.1 and the prediction for the 
rapidity dependence of the Rd(z) for three rapid¬ 
ity regions chosen to match the rapidity regions 
used in the measurement of the jet Raa by AT¬ 
LAS [11], namely |p| < 0.3, 0.3 < \y\ < 0.8, and 
1.2 < IpI < 2.1. The differences in the Rd{z) are 
better quantified in terms of the ratio of Rd{z) eval¬ 
uated in a given rapidity region to that evaluated in 
the region of |p| < 2.1 as shown in the right panel of 
Fig. 11. The ratio is predicted to reach larger val¬ 
ues in the more central rapidity region and smaller 
values in more forward region compared to the in¬ 
clusive rapidity interval. Both of these effects are in 
the maximum at the level of 6-7% of the inclusive 
Rd{z)- This predicted behavior represents another 
possibility to test the model and, more generally, 
to probe the differences between the energy loss of 
quark and gluon jets. 

7. Modeling Z3(px) and charged particle 
Raa 

We have shown that our simple model for the 
medium modifications of single jets can reproduce 
the measured nuclear modification factor and its ra¬ 
pidity and transverse momentum dependence. The 
model can also reproduce some of the qualitative 
features seen in measured inclusive jet fragmenta¬ 
tion functions, namely the suppression at interme¬ 
diate z and an enhancement at large z. It was 
shown that these features in the measured fragmen¬ 
tation functions arise from the change in the quark 
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Figure 9: Ratios of D(z) distributions for six bins in collision centrality to those in peripheral (60-80%) collisions, 

measured by ATLAS for R = 0.4 jets [12] (black markers) are compared to the analytic calculation 
(red line) and MC calculation (blue histogram) of the same quantity in the non-constant fractional energy loss model. 


(or gluon) fraction due to the medium-induced en¬ 
ergy loss. Since the model can explain both the jet 
spectra and the fragmentation functions, it should 
be also able to reproduce the charged particle trans¬ 


verse momentum distribution, of charged 

particles produced within jets, and the nuclear 
modification factor of charged particles, i?AA) mea¬ 
sured at high transverse momenta of charged par- 




Figure 10: Predicted jet Rw as a function of pj);* in forward 
rapidity intervals 2.1 < \y\ < 2.8 (left) and 2.8 < \y\ < 3.5 
(right) shown for four centrality bins, 0-10% (black), 20- 
30% (blue), 30-40% (red), and 60-80% (green). See text for 
details. 
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Figure 11: Left panel: ATLAS data on the ratio of D{z) distributions, Rd{z)-> 0-10% to 60-80% centrality bin (black 

markers) compared to the analytic calculation of the non-constant fractional energy loss for three different values of color 
factors: the default value of color factor Cp = 9/4 (blue), Cp = 1.0 (red), and Cp = 3.0 (green). Middle panel: Prediction of 
the ratio of D(z) distributions, Rd{^z)-> 0-10% to 60-80% centrality bin in different jet rapidities. The Rp){^z) the rapidity 

region of \y\ < 2.1 (black) is compared to the Rp)(^z) ^^le regions: \y\ < 0.3 (green), 0.3 < \y\ < 0.8 (blue), and 1.2 < \y\ < 2.1 
(red). Right panel: The ratio of predicted Rd(^z) different rapidity regions to the Rp>(^z) region of \y\ < 2.1. 


tides, p^. The D{p^) distributions were mea¬ 
sured by ATLAS [12] and CMS [13] for jets with 
> 100 GeV. The charged particle at 
high p^ was measured by CMS [8]. Compared to 
the inclusive jet fragmentation functions which are 
largely independent of the jet transverse momen¬ 
tum for a given flavor of the initial parton, these 
observables couple together the change in the jet 
fragmentation and the change in the underlying jet 
spectra. Thus, these observables provide another 
important input to test the model. 

The PYTHIAS simulated events were used to 
simulate the modifications of the D{p^) distribu¬ 
tions. The simulation was done in two steps. In 
the first step, the D{p‘^) distributions were booked 
in I GeV bins of the p!^^ to two look-up tables, for 
quark initiated and gluon initiated jets separately. 
In the second step, the jet suppression was applied 
at the single jet level in the same way as for the case 
of simulating the inclusive jet Raa- For each sup¬ 
pressed jet of a given flavor, the D{p^) distribution 
corresponding to the quenched jet px was read off 
from the look-up table and added to the histogram 
which formed the resulting D{p^) distribution af¬ 
ter being properly normalized by a total number 
of quenched jets with > 100 GeV. The central 
to peripheral ratio of D{p'i^) distributions, R^x^p^p 


was evaluated. The result is shown along with the 
data by ATLAS in Fig. 12. The figure shows that 
our model can reproduce the qualitative features 
observed in the data, namely the suppression of 
yields at intermediate plj? and an enhancement at 
high pljk. This is not surprising given the success 
of the model in describing the D{z) modifications, 
but it represents an important consistency check. 

To evaluate the nuclear modification factor of 
charged particles, R^aa^ same pro¬ 

cedure as for evaluating the Ru(p^) with only two 
differences which is evaluating the D{p‘^) distribu¬ 
tions for jets with no threshold on jet px and avoid¬ 
ing a normalization of the final D{p'^) distributions 
by the total number of jets. This approach is based 
on the fact that each charged particle with a given 
plj? must come from the jet with pl^* > plj?. This al¬ 
lows to construct the R'aa p't ^ GeV which 
is the kinematic region where we have a good con¬ 
fidence in modeling the inclusive jet suppression. 
The resulting R'aa shown along with the R'aa 
measured by GMS in Fig. 13. 

The figure shows that the model can reproduce 
the qualitative features seen in the data at high- 
p^, namely the increase of the R'aa with increas¬ 
ing pJjk and its centrality dependence. However, the 
slope of the charged particle Raa in the model dif- 
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fers from that in the data at the lowest values 
included in this analysis and the model systemat¬ 
ically slightly over-predicts the i?AA the data. 
These disagreements may be due to insufficient pre¬ 
cision in the modeling of distributions by 

PYTHIAS, something that can and will be tested in 
future analyses. There may be consequences from 
the disagreements between the data and the model 
in the Rd{z) including the low -2 (low- 

p'^) excess and fact that the data is systematically 
lower than the model in the range 0.1 < z < 0.2 
(10 < < 20 GeV). The former is discussed in 

the next section, the latter may be insufficient to 
explain the observed difference. Generally, the dis¬ 
agreement may be due to jet quenching physics not 
included in the model. In particular, the possibility 
that the jet quenching produces changes in the frag¬ 
mentation functions of lower-pl^* jets different from 
that observed in the > 100 GeV jets cannot be 
excluded. 

8. D{z) low-z excess 

As described above, our model for the modi¬ 
fication of the fragmentation function in Pb-|-Pb 
collisions is based on the assumption that energy 
lost during the evolution of a parton shower in the 
medium does not appear as part of the measured 
jet and that the final fragmentation products have 
the same D{z) distribution as if the reduced-energy 
jet fragmented in vacuum. While the validity of 
these assumptions may be debated (next section), 
the soft excess in the Pb-fPb D{z) distributions, 
which cannot be explained by the different quark 
and gluon quenching, presents a manifest violation 
of the assumptions of the model. Supposing that 
the enhancement at low-z reflects either radiation 
within the jet, recoil partons, or collective response 
of the medium, from the point of view of our model, 
the extra low-z hadrons provide an upward shift of 
the jet energy. In the ATLAS Pb-l-Pb fragmenta¬ 
tion function measurement, the excess fragments in 
the range 0.02 < z < 0.04 were found to contribute 
~ 2% of the jets transverse momentum for the 0- 
10% centrality bin. Since that estimate could not 
account for contributions from hadrons below the 
minimum transverse momentum of the charged par¬ 
ticle measurement, it is clearly an under-estimate 
of the contribution of “excess” low-z partons to the 
jets energy. However, we take this number to be an 
order-of-magnitude estimate of the fractional con¬ 


tribution of the low-z excess to the energy of the 
typical jet which we will refer to as 

Since the energy in the low-z fragments likely is 
proportional to the jets energy, we assume that 

^soft ^ /int^soft ^ _ /int)$|oft^ (20) 

with and where <1)®°^*, and 

represent the average soft excess in quark, 
gluon, and all jets, respectively. The contribution of 
the excess soft hadrons to the jet energy is, in prin¬ 
ciple, already effectively accounted for in our anal¬ 
ysis of the jet suppression where it will reduce S. 
However, it could affect our description of the mod¬ 
ified fragmentation functions by increasing the jet 
energy that appears in the denominator of the z def¬ 
inition in the data by a fraction l-|-<i)®°^*, and, thus, 
reducing the z values by a factor 1/(l-f 
Thus, the modification of the fragmentation func¬ 
tion in the data at z should correspond to the mod¬ 
ification in our model at a z value of z (l -I- 
To reproduce that effect, we have corrected Eq. 14 
as follows, 

^meas(^) ^ ^ + 

(l-/“‘)i^,(z[l + cl>ft])^ 

and recalculated the D(z) modifications taking 
^inc* 0.02. The results are shown in the left 
panel of Fig. 14. Accounting for the shift in the 
jet energy shifts Rd{z) down and leads to a better 
agreement between the data and the model. More 
than one sigma disagreement my now be seen only 
for the data point near z = 0.1. 

The depletion seen in the modeled Rd{z) at high- 
z reflects the rapid decrease in the parameterized 
D(z) for z > 1. The result in Fig. 14 is continuous 
across z = 1 for reasons given in Sec. 2, but the 
D{z) still falls rapidly near z = 1. The data do not 
yet have the precision to resolve a change in behav¬ 
ior near z = 1 like that shown in the figure. Future 
measurements with improved precision could test 
for such an effect. Observing the depletion would 
provide strong empirical support to the picture that 
(most) jets fragment as in the vacuum but with ad¬ 
ditional energy from low-z particles whose origin is 
not yet understood. 

The contribution of the excess soft hadrons to the 
jet momentum also influences the modeled Ro{p^) 
and The impact of this contribution can be 

avoided by reducing the quenched jet momentum 
used to look-up for the D{p^) distributions corre¬ 
sponding to quenched jets by The impact of 
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Figure 12: Ratios of D{p^) distributions for six bins in collision centrality to those in peripheral (60-80%) collisions, 
D(p^)|cent/^(p^)|60—80? measured by ATLAS for R = 0.4 jets [12] (black points) are compared to the MC calculation 
(blue histogram) of the same quantity in the non-constant fractional energy loss model. 
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Figure 13: Nuclear modification factor of charged particles measured by CMS [8] in six centrality bins (black points) is 
compared to the MC calculation of the same quantity in the non-constant fractional energy loss model (blue histogram). 
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Figure 14: Demonstration of the effects of accounting for the contribution of low- 2 : particles in jet energies (see text). Left: 
calculated R]j{^z) with (red) and without (blue) the low -.2 excess; middle: same for -Rd(px)’ with and without the 

soft- 2 : excess. 


this on modeled and = 0.02 is 

shown in the middle and right panel of Fig. 14, re¬ 
spectively. While the impact of soft excess hadrons 
on the distribution is rather small the im¬ 

pact on R^^p^ is more significant leading to a better 
agreement of modeled i?AA with the data. 

9. Discussion 

This paper has presented an analysis of impli¬ 
cations of recent data on single jet suppression 
and inclusive jet fragmentation in Pb+Pb colli¬ 
sions. The analysis was based on a simple model for 
the quenching of a parton shower in the medium, 
namely that the parton shower loses energy to the 
medium in a manner such that the lost energy does 
not appear within the jet “cone”. The fragmenta¬ 
tion of the resulting jet was assumed to be the same 
as the fragmentation of a jet in vacuum. These as¬ 
sumptions may seem unreasonably simplistic given 
the current understanding of the evolution of high- 
energy parton showers in vacuum where pertur¬ 
bative calculations can describe many features of 
the resulting distributions of fragments. However, 
studies of the impact of the medium on the color- 
coherence of parton showers [25] suggest that the 
medium is unable to resolve the internal structure 
of many jets. According to Ref. [26] those jets in¬ 
teract with the medium as if they consist of a single 
color charge, and the reduced-energy jet fragments 
as if it were in vacuum. An extensive analysis of the 


angular and longitudinal momentum distribution 
of the radiated energy [26-31] indicates that the 
medium-induced radiation flow to large angles con¬ 
sistent with CMS measurements [32] . The combina¬ 
tion of these two idea/results suggests a picture for 
jet-medium interactions that is qualitatively similar 
to that used in this analysis. 

This work started from a simple hypothesis that 
the dependence of the measured jet i?AA and 
the modifications to the jet fragmentation functions 
could both be explained by the different quenching 
of quarks and gluons and the dependence of the 
primordial quark fraction. Indeed, the first result 
using analytic expressions, power-law spectra and a 
constant fractional shift provided remarkably good 
agreement with both the Raa and the D(z) data. 
However, that success was short-lived as it was 
found to disagree with a Monte-Carlo implementa¬ 
tion that used the PYTHIAS quark and gluon spec¬ 
tra directly. When deviations of the spectrum from 
the pure power-law form were accounted for, the 
discrepancy between analytic and Monte-Carlo re¬ 
sults was resolved, but the constant-fractional shift 
could no longer describe the data. The sensitivity of 
the jet Raa to the shape of the jet spectrum is well 
known, but this analysis provides clear demonstra¬ 
tion of the sensitivity of interpretations of the Raa 
measurements to the accuracy in the description of 
the primordial jet spectra. 

The analysis presented in this paper appears to 
rule out the possibility of a constant fractional shift 
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parameterization of the effects of quenching on the 
jet spectrum. While a constant fractional shift is in¬ 
compatible with most energy loss calculations, the 
weak dependence of the measured jet suppres¬ 
sion has often been (informally) interpreted as indi¬ 
cating fractional energy loss. However, as shown in 
Fig. 4, even when accounting for the fact that /q in¬ 
creases with a constant fractional energy loss 
produces an Raa that decreases with due to 
the fact the the primordial jet spectrum steepens 
relative to a pure power law with increasing pl^*. 
The pl^*-dependent shift extracted from the data, 

in fact, varies like as predicted in Ref. [21]. 

That reference makes clear that the shift should not 
necessarily be interpreted as the average energy loss 
of the jets. Indeed, the shift approximately param¬ 
eterizes the convolution of the energy loss distribu¬ 
tion with the primordial jet spectrum. Thus, direct 
interpretation of the s' values extracted in this anal¬ 
ysis is not immediately possible. Nonetheless, with 
the simple centrality dependence of s' observed in 
Fig. 7, the full pJ^* and centrality dependence of the 
jet suppression can be accounted for by three pa¬ 
rameters, a and the effective slope and intercept of 
the A^part dependence of s'. 

The apparent accidental cancellation of the vari¬ 
ation of spectrum shapes and quark fraction ob¬ 
served here and in Ref. [33] that leads to the mea¬ 
sured constancy of Raa with rapidity [ 11 ] is unfor¬ 
tunate as it reduces the sensitivity of the measure¬ 
ment to features of the energy loss that could help 
to test or constrain calculations of medium-induced 
energy loss. However, the extensions of our analysis 
to larger rapidity indicate that at larger rapidities, 
the increased curvature of the primordial jet spec¬ 
tra will provide a different and stronger variation 
of Raa with pl^*. That is particularly true in the 
2.8 < IpI < 3.5 bin where measurements may be 
difficult in the most central collisions but could be 
performed in more peripheral centrality bins where 
the effect should still be measurable. 

The analysis in this paper suggests that much or 
all of the observed modifications to the jet fragmen¬ 
tation functions in Pb-f Pb collisions can arise from 
the different quenching of quark and gluon jets - ex¬ 
cept ^or the enhancement at low z. This observation 
could have important implications for the theoret¬ 
ical understanding of the quenching physics. For 
example, in the context of the picture presented in 
Ref. [25] where quenching of the jets is influenced 
by the ability of the medium to resolve the internal 


structure of the jet, jets that are not resolved by the 
medium fragment according to their vacuum frag¬ 
mentation functions. But, the measured D{z) dis¬ 
tributions will still differ from those in pp or periph¬ 
eral Pb-|-Pb collisions due to the different energy 
loss of the quarks and gluons. On the other hand, if 
the fragmentation functions of the quark and gluon 
jets are separately modified, then we would have to 
conclude that the medium is resolving the internal 
structure of the jet. The simplicity of the analysis 
used in this paper is not adequate to draw a firm 
conclusion that the D[z) modifications at interme¬ 
diate and large z can be completely attributed to 
flavor-dependent quenching. However, the analy¬ 
sis shows the importance of explicitly addressing 
the effects of different quark and gluon quenching 
in future analyses of the fragmentation functions. 
Measurements of jet fragmentation in 7 -jet events, 
where the jet spectrum has a larger quark fraction 
would be valuable in addressing this issue. 

The analysis in this paper suggests that the 
magnitude and transverse momentum dependence 
of charged particle Raa for particles with px > 
20 GeV can largely result from the different quench¬ 
ing of quark and gluon jets as well. The disagree¬ 
ment between the measured charged particle Raa 
and the Raa in the model, namely larger suppres¬ 
sion and steeper Raa as a function of px seen in the 
data, remains to be understood. It may be arising 
from missing physics in the model, difference be¬ 
tween the measured charged particle spectra and 
their PYTHIAS simulation, or from larger contri¬ 
bution of soft particles to quenched jets at lower 
Px. Irrespective of the source of this disagreement, 
the result clearly points to the importance of un¬ 
derstanding the measurements involving single par¬ 
ticles in the context of fully reconstructed jets. 

As described above, the enhanced production of 
hadrons at low z observed in data [12, 13] cannot 
be explained by the different energy loss of quarks 
and gluons. It’s also interesting that it could not 
be explained by a strong-coupling calculation of en¬ 
ergy loss [34] but did arise in a collisional energy- 
loss scenario also tested in that same paper, sug¬ 
gesting that the excess could arise from recoiling 
constituents of the medium [35]. If this explana¬ 
tion is correct, then the low-z excess is directly 
probing (part of) the medium-response to the pas¬ 
sage of jets. An alternative explanation was pro¬ 
vided in Ref. [31] that showed a soft-z excess can 
arise when the medium resolves the constituents of 
the jet core. More speculatively, it might be that 
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the excess arises from the collective response of the 
medium such as a diffusion wake [36]. Regardless 
of the explanation, the energy contributed to the 
low-z particles acts in the context of our model 
as an extra contribution of the energy of the jet. 
The effect of this contribution produces modest but 
noticeable effects on the fragmentation functions. 
We note that the persistence of the low-z excess in 
non-central collisions indicates that it does not re¬ 
sult from systematics in the measurement. Given 
the potential importance of the low- 2 ; excess, more 
detailed measurements in non-central events where 
the effects of the underlying event are smaller may 
be warranted. 

10. Conclusions 

This paper has presented an analysis of single jet 
measurements from Pb-|-Pb collisions at the LHC 
using a simple, phenomenological model. The anal¬ 
ysis used quark and gluon spectra obtained from 
PYTHIAS and applied the “shift” formulation de¬ 
scribed by Baier et al. to describe the single-jet 
suppression. The modifications to the jet fragmen¬ 
tation functions were assumed to result from the 
different quenching of quarks and gluons assuming 
that the quenched jets fragment as they do in vac¬ 
uum. 

Our analysis showed that the transverse momen¬ 
tum dependence of the quark fraction plays a role 
in the evolution of the jet Raa with as was 
also observed in Ref. [33]. However, the curvature 
of the primordial jet spectrum relative to a pure 
power-law also substantially effects the Raa- The 
data were found to be incompatible with a constant 
fractional shift, S = Fits of the data using 

a shift that varied with as yielded a 

value a ~ 0.55 compatible with the prediction in 
Ref. [21]. The measured modifications of the inclu¬ 
sive jet fragmentation functions can be largely ex¬ 
plained as resulting from the different quenching of 
quarks and gluons. More specifically, the suppres¬ 
sion of the fragmentation function at intermediate 
z and the enhancement in the fragmentation func¬ 
tion at large z reflect the different shapes of the 
quark and gluon fragmentation functions and the 
increase in the quark fraction due to the greater 
energy loss of gluon jets. However, our model is 
not able to account for the low-z enhancement in 
the Pb-|-Pb fragmentation functions. If our analysis 
is correct, the fragmentation function modihcations 


provide a direct test of the color-charge dependence 
of jet quenching. The model can explain most of the 
observed suppression of the charged particle yield 
at high plj!' though it slightly under-predicts the 
amount of suppression and has a weaker de¬ 
pendence at the lowest values included in the 
analysis. 

Because the quark fraction varies as a function 
of rapidity, both the Raa and the D{z) modifica¬ 
tions should evolve as a function of rapidity, though 
the rapidity-dependent evolution of the quark and 
gluon spectra appears to cancel out the effects of 
the changing quark fraction over the rapidity range 
measured in the ATLAS data. Predictions are 
made for Raa and Rd(z) at large rapidity that 
would allow the conclusions of this analysis to be 
tested. 
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